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ABSTRACT 


New  developments  on  lifting  surface  theory  for  oscilla¬ 
tory  subsonic  flow  are  given.  Key  concepts  in  the  analysis 
are  the  use  of  concentrated  loads  rather  than  distributed 
pressure  "mode"  shapes,  and  the  development  of  a  modified 
kernel  function  which  gives  average  values  of  vertical  velo¬ 
city  over  chosen  intervals  and  which  eliminates  all  singular¬ 
ity  problems. 

Features  of  the  analysis  are:  (l)  loads  are  given  directly 
in  terms  of  vertical  velocities,  (2)  no  pressure  modes  have  to 
be  assumed,  (3)  singularities  are  obviated,  (4)  the  locations 
of  control  downwash  points  are  specified  systematically, 

(5)  control  surfaces  may  be  included,  (6)  treatment  of  non- 
planar  surfaces,  such  as  P-tails,  is  possible,  and  (7)  appli¬ 
cation  is  made  through  a  simple  quick  routine  procedure. 
Examples  are  given  throughout  to  illustrate  the  concepts. 


This  abstract  is  subject  to  special  export  controls,  and 
each  transmittal  to  foreign  governments  or  foreign  nationals 
may  be  made  only  with  prior  approval  of  the  Air  Force  Flight 
Dynamics  Laboratory  (FDTR),  Wright-Patterson  Air  Force  Base, 
Ohio  45433. 
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speed  of  sound 
see  equation  (66a) 
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velocity  components  in  the  z -direction 
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coordinates  to  center  of  an  area  element 
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angle  of  attack;  also  a  =  e 
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interval  in  the  x-directlon 
alternate  coordinate  in  the  y-direction 
interval  in  the  y-direction 

alternate  coordinate  in  the  x-dlrectlon 
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velocity  potential 

strength  of  potential  sheet 
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surfaces  of  potential  sheet 
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Notations: 
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1(x  “  xr)  >  unit  step  function  at 


x  =  x 
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,  Dirac  function  at  y  =  0 


SECTION  I 


INTRODUCTION 


This  report  Is  concerned  with  oscillatory  airforces  but 
is  motivated  by  the  problem  of  determining  the  frequency 
response  function  of  an  aircraft.  The  frequency  response 
function,  defined  in  a  general  sense  as  the  structural  response 
due  to  a  unit  sinusoidal  input,  is  now  an  Important  ingredient 
in  the  treatment  of  some  of  the  major  aeroelascic  problems  of 
an  aircraft.  Design  for  gusts  by  power  spectra^  techniques  is 
one  area.  Flutter,  representing  a  special  situation  of  the 
frequency  response  function,  is  another.  Because  of  this  prom¬ 
inence,  a  desire  exists  to  establish  the  function  in  a  quick 
reliable  way. 

The  establishment  of  the  frequency  response  function  is 
conveniently  broken  down  into  two  areas,  one  dealing  with  the 
aircraft  structure,  the  other  with  the  nonsteady  aerodynamics. 

The  structure  is  generally  handled  by  either  of  two  approaches, 
one  involving  the  use  of  modal  functions,  the  other  the  use  of 
lumped  or  discrete  masses.  The  aerodynamics  is  usually  handled 
by  modified  strip  theory  or  by  a  lifting  surface  approach 
involving  the  use  of  a  kernel  function. 

In  appraising  the  problem  of  determining  the  frequency 
response  function,  it  appeared  that  'jnpr ovement s  in  the  handling 
of  the  aerodynamics  were  needed.  These  possible  improvements 
are  discussed  here  mainly  with  reference  to  the  use  of  lifting 
surface  theory.  In  spite  of  the  fact  that  this  theory  has  been 
applied  successfully  in  various  instances,  some  shortcomings 
still  exist.  The  approach  does  not  seem  to  be  as  computation¬ 
ally  economic  as  is  strip  theory  (because  of  this  fact,  modified 
strip  is  still  used  in  the  treatment  of  very  large  flexible  air¬ 
craft).  Pressure  "modes"  have  to  be  assumed  from  which  the 
downwash  is  computed;  a  desire  is  to  be  able  to  determine  the 
pressure  directly  from  a  prescribed  downwash  condition.  Uncer¬ 
tainty  still  exists  on  how  many  collocation  points  should  be 
used,  and  where  they  should  be  located.  The  technique  for 
including  control  surfaces  is  not  clear;  only  recently  have  some 
special  studies  of  control  surfaces  been  made.  Finally,  lifting 
surface  theory  has  been  restricted  mainly  to  planar  surfaces. 

The  purpose  of  this  report  is  to  develop  improvements  in 
the  use  of  lifting  surface  theory;  and,  in  turn,  to  ease  the 
problem  of  determining  the  frequency  response  function.  Its 
aims  are  to:  (l)  allow  loads  to  be  established  directly  in  terms 
of  given  downwash  values,  (2)  eliminate  the  necessity  of  assuming 
pressure  modes,  (3)  make  the  choice  cf  downwash  points  less  arbi¬ 
trary,  (4)  include  control  surfaces,  (5)  include  nonplanar  sur¬ 
faces,  and  (6)  simplify  the  application.  Key  concepts  in  the 
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development  are  the  use  of  concentrated  loads  rather  than  dis¬ 
tributed  pressure  mode  shapes  and  the  development  of  a  modified 
kernel  function  which  gives  average  values  of  vertical  velocity 
over  chosen  intervals  and  which  eliminates  the  problem  of  deal¬ 
ing  with  singularities. 

The  present  treatment  is  confined  to  the  subsonic  range, 
but  many  of  the  ideas  should  apply  equally  well  to  the  super¬ 
sonic  range  and  especially  in  the  treatment  of  problems 
involving  mixed  subsonic  and  supersonic  flow. 


SECTION  II 

GOVERNING  DIFFERENTIAL  EQUATIONS  AND  BOUNDARY  CONDITIONS 


Governing  Differential  Equations.-  The  developments  given 
in  this  report  could  start  directly  with  the  well  established  and 
generally  accepted  linear  equations  for  3-dimensional  oscillatory 
subsonic  flow.  For  completeness,  however,  a  brief  review  of  their 
derivations  is  given  in  this  section. 


We  consider  a  lifting  surface  immersed  in  a  flow  from  the 
left  as  shown  in  the  following  sketch. 


Wi 


The  flow  is  considered  to  be  isentropic. 
differential  equations  are: 

Euler  eauations. 


The  governing 


Equation  of  continuit 


PlVl  +  PlWl 


Isentroplc  gas  law  (from  equation  of  state), 

Pi 


=  c 


(3) 


Pi 


We  consider  a  small  perturbation  state  by  writing 


U1 

=  U  +  u 

V1 

=  V 

W1 

=  w 

Pl 

Q. 

+ 

o 

Q. 

II 

Pl 

=  Po  +  p 

where  u,v,w,p,  and  p  are  the  perturbations  which  are  small 
relative  to  the  steady  state  terms.  The  substitution  of 
equations  (4)  into  equations  (1),  (2)  and  (3)  yields,  after 
discarding  the  second  order  terms,  the  following  linearized 
equations 


3T  +  u 

37 

P0  5^ 

u 

d 

37 

i  a 
p0 

^  1  TT 

a 

1  d 

+  U 

37 

a 

3 

1 

I'd  ,  „ 

37 

37 

“2— 
a  p„ 

+  U  c 

u 


w 


=  0 


(5) 


3x/ 


The  speed  of  sound  a  appears  in  the  last  equation  as  a  result 
of  using  equation  (3)  to  replace  p  in  equation  (2)  by  p  ; 
specifically  equation  (3)  indicates 


dpl  £  pl 

a?T  p  = 


(6) 
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The  treatment  of  equations  (5)  is  simplified  by  introducing 
the  velocity  potential  0  such  that 


With  these  equations  all  of  the  first  three  equations  of  equations 
(5)  reduce  identically  to  the  equation  for  pressure 

?  “  '  p0(u  H  +  M)  (8) 

while  the  last  equation  becomes  the  governing  differential 
equation 


Equations  (8)  and  (9)  are  the  basic  well  known  equations  that 
are  to  be  treated  in  this  report .  It  is  significant  to  note 
that  equations  (5)  show  that  equation  (9)  also  applies  when 
any  of  the  independent  variables  u,v,w,  or  p  is  written  in 
place  of  0  . 

Boundary  conditions.-  The  problem  to  be  solved  is  fixed 
by  stating  the  boundary" conditions  that  must  be  used  in  con¬ 
junction  with  equation  (9).  In  developing  the  boundary  con¬ 
ditions  It  is  instructive  first  to  examine  the  symmetry  pro¬ 
perties  of  the  variables  u,v,w,p  and  0  .  We  note  that  the 
vertical  velocities  on  the  upper  and  lower  surface  of  the  wing 
must  be  equal  (wing  considered  thin)  and  therefore  reason  that 
w  must  be  a  symmetrical  function  with  respect  to  z  ;  equations 
(5)  and  (8)  indicate  then,  that  u,v,p  and  0  must  all  be  anti- 
symmetrical  with  respect  to  z  .  These  symmetry  properties  allow 
the  flow  problem  to  be  solved  by  considering  the  upper  half 
z-region  only.  Boundary  conditions  on  the  x  -  y  plane  thus 
become  of  concern.  It  is  convenient  to  specify  these  boundary 
conditions  in  terms  of  the  three  different  regions  A,  S,  and  W 
shown  in  figure  1.  In  the  region  A  ahead  of  the  wing,  the 
only  perturbations  that  can  exist  are  those  propagated  forward; 
no  discontinuities  can  arise,  and  thus,  because  0  is  anti¬ 
symmetric,  0=0  on  the  x  -  y  plane.  For  the  region  S  , 
the  wing  supports  a  pressure  discontinuity  across  its  thickness 
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with  p  =  -  p~  ,  where  +  and  -  denote  respectively  the  upper 
and  xower  surfaces.  The  w  in  this  region  must  be  that  given 
by  the  wing.  The  region  W  is  established  by  the  Kutta  condi¬ 
tion  whit h  states  that  w  must  remain  finite  at  the  trailing 
edge,  or  alternatively  that  p  =  0  at  this  edge.  From  tnis 
condition  it  is  reasoned  that  the  wing  leaves  behind  it  a  dis¬ 
continuous  wake  sheet,  which  must  be  continuous  in  w  ,  which 
cannot  support  a  pressure  jump,  and  which  therefore  must  be 
discontinuous  in  u  ,  v  ,  and  0  .  Figure  1  specifies  the 
problem  completely.  The  main  difficulty  is  the  problem  of 
dealing  with  mixed  boundary  conditions. 

Throughout  the  remainder  of  the  report  we  will  be  concerned 
with  the  oscillatory  case  only,  where  the  independent  variables 
0  and  p  may  be  expressed  in  the  form 


0 


0(x,y,z)e 


iajt 


P 


p(x,y,z)e 


iajt 


and  similarly  for  u,  v,  w  and  p  .  Thus  for  the  oscillatory 
case,  equations  (9)  and  (8)  become 


(10) 


p  =  -  p0(u  If  +  ic#) 

The  harmonic  term  eiaJt  will  be  suppressed  for  brevity  in 
writing  in  all  equations  to  follow.  As  mentioned  earlier, 
equation  (10)  applies  not  only  for  0  but  for  the  other 
independent  variables  as  well;  thus 

V2p  — ^  +  ico)  p  =  0  (12) 

E 


The  integration  of  equation  (11)  gives  a  result  which  is 
basic  to  all  subsequent  development;  this  integration  yields 

ictac  x 

0(x,y,z)  =-ie  ~U~  I  e  (13) 

Ho  u 

—oo 
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Of  concern  also  is  the  result  obtained  from  this  equation  as 
z  -*  0  ;  this  result  is 


iaac  x 

0+(x,y,o)  =  -  -pjj  e  U  Te  U  p+(£,y,o)d£  (14) 

^o  u 

-00 


where  the  +  sign  signifies  values  at  the  upper  surface  of  the 
discontinuous  sheets.  Equations  (10)  through  (14)  form  the 
basis  for  the  developments  in  the  subsequent  sections.  In  a 
general  sense,  solutions  are  made  by  means  of  equations  (13)  or 
(1 4) ,  in  terms  of  basic  source  solutions  of  either  equations 
(10)  or  (12),  and  with  due  regard  being  given  to  the  boundary 
conditions . 

Basic  source  solutions  of  equations  (10)  arid  (12).-  Some 
solutions  of~equations  (10)  and  (12)  that  form  the  base  for  the 
solutions  of  various  oscillatory  flow  cases  are  the  following 
(see  reference  2): 

Monopole  (of  unit  strength). 


-ico 

0m  or  pm  =  ‘  TF®  6  aP 
Dipole  (of  unit  strength), 

p  -io>  ^*1* 

OT  Pd  =  '  t^3  (X  +  ^ 

where 

R  =v^2  +  p2(yS  +  z2) 


(15) 


(16) 


These  two  equations  are  related  by  the  operation 

d0_ 

j  _  m 

*a  ~  "  H" 
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Two  other  solutions  that  are  used  herein  are  the  following 
reduced  cases: 


Line  dipole  for  2 -dimensional  Incompressible  flow. 


0  or  p 


1  z 


(17) 


Dipole  for  3-dlmensional  incompressible  flow. 


0  or  p 


1  z 


R 


.  c_  .  c 

+  y  +  z 


(18) 


Point  function  characteristics  of  Eqs.  (16),  (17),  (18) . - 
We  remark  here  on  the  point  function  character  of  equations  (16) 
(17)  and  (18),  especially  with  respect  to  the  use  of  a  factor  2; 
some  concern  arises  at  times  as  to  whether  or  not  a  factor  2 
should  be  included,  and  these  remarks  are  offered  for  possible 
clarification.  Consider  a  2-dimensional  pressure  discontinuity- 
sheet  extending  from  x  =  x^  to  x  =  X2  .  The  strength  is 
specified  by  the  function  pQ(x)  ,  which  by  the  antisymmetric 
nature  of  p  implies  a  pressure  of  -  l/2p0  on  the  upper  sur¬ 
face  of  the  sheet  and  +  l/2p0  on  the  lower  surface.  The  di¬ 
pole  strength  in  an  elemental  length  dx  is  p0(x)dx  ,  and  thus 
from  equation  (17)  the  field  pressure  may  be  found  as 


p(x,z) 


zP0U)d£ 
(x  -  £)*  v  z‘ 


(19) 


In  considering  this  equation  in  the  limit  z  -*■  0  ,  we  note  that 
only  the  portion  of  the  integrand  in  the  immediate  vicinity  of 
x  can  contribute  to  the  integral.  The  equation  may  be  reduced 
therefore  to 


p(x,z) 


ZdTJ 


+  z 


2 


P0(£)  may  be  considered  constant. 


This  equation  yields 


p+(x,o)  =  11m  — ~  —  tan"1  D.  I 

‘-IT  Z  1 

z  -*  o  -e 


2  P0(x)  (20) 

thus  indicating  the  pressure  to  be  the  upper  surface  pressure 
?hus  alsoUsei„?he  P°lnt  funotlon  character  of  equation  (17)  Is 

equation  .fSS.**  SqUatl0n  (20)  15  subs«tuted  in 


p(x,z) 


-hi 


*  ( y  4\2  j.  2 

(X  “  *)  +  Z 


(21) 


and  we  thus  make  the  following  rule* 

2  dlacontlnuity  strength  is  specified  use  equationS(17) 
as  is  to  deduce  the  field  pressure  (equation  (IQ))  •  if  the  umv»r 
surface  pressure  is  specified  use  equation  (17)  with  the  facto? 
suppressed  to  deduce  the  field  pressure  (equation  21)) 
observations  may  seem  trivial,  but  trouble  has  occuriid 
in  instances  where  they  have  not  been  fully  understood. 

h  •,  j  f  similar  point  function  proof  and  rule  may  be  shown  to 
hold  in  the  application  of  equations  (16)  and  (18). 
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SECTION  III 


3QLUTION  BY  NUMERICAL  ORIDWORK  SCHEMES 


The  oscillatory  problem  under  consideration  here  may  be 
stated  as  follows.  Consider  a  wing  to  be  oscillating  sinusoid¬ 
ally  in  the  z-direction  in  some  prescribed  deflection  shape;  we 
wish  to  find  the  particular  solution  of  equations  (10)  or  (12) 
which  gives  the  pressure  distribution  on  the  wing  and  for  which 
all  the  essential  boundary  conditions  are  satisfied.  Such  a 
direct  so\ution  has  been  possible  for  only  a  few  special  cases, 
such  as  2-dimensional  incompressible  flow.  For  the  more  general 
cases  solution  has  had  to  proceed  in  an  inverse  sense;  thus,  a 
pressure  distribution  is  assumed  and  solution  is  made  for  the 
vertical  velocities  that  are  induced.  To  be  more  specific  let 
us  review  the  basic  steps  of  the  kernel -function  approach  that 
has  been  used  in  the  past.  Equation  (16)  for  p^  is  sub¬ 
stituted  into  equation  (13) »  from  which  the  vertical  velocities 
are  found  by  equation  (7c);  formally  the  result  appears  as 


w  =  K(x,y,o)  = 
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This  function  is  called  the  kernel  function;  it  is  analogous 
to  an  influence  coefficient  development  and  specifically  denotes 
the  vertical  velocity  that  is  developed  at  the  point  (x,y,o) 
due  to  a  unit  pressure  dipole  at  the  origin.  Its  evaluation 
and  a  discussion  of  the  associated  singularities  are  given  in 
reference  1.  With  this  equation  the  vertical  velocity  that 
results  from  a  given  loading  distribution  pQ  over  the  wing 
follows  as 


w 


=  ff  P0(M)  K  (x  -  e,y 
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(23) 


To  insure  satisfaction  of  the  Kutta  condition,  the  distribution 
is  chosen  so  that  pn  is  zero  at  the  trailing  edge.  To  obtain 
the  solution  for  specific  lifting  surfaces  the  technique 
generally  used  is  to  introduce  several  pressure  "modes"  or 
distributions  and  then  by  a  collocation  technique  to  find  what 
combination  of  these  pressure  modes  leads  to  the  desired  vertical 
velocities  on  the  wing.  In  terms  of  the  wing  deflection  Z  ,  and 
in  terms  of  any  gust  velocities  Wg  that  also  may  be  under  con¬ 
sideration,  these  vertical  velocities  are  given  by 
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where  the  second  relation  is  for  the  sinusoidal  case.  It  is 
noted  that  the  evaluation  of  K  and  its  use  in  equation  (23) 
is  fraught  with  horrible  singularities. 


The  treatment  to  follow  has  three  objectives:  (1)  to  make 
a  direct  solution  of  the  pressures  in  terms  of  the  vertical 
velocities,  (2)  to  obviate  the  necessity  for  introducing 
pressure  "mode  shapes",  and  (3)  to  greatly  simplify  or  eliminate 
the  problem  of  dealing  with  the  difficult  singularities.  Four 
possible  procedures  are  advanced. 

Scheme  I. -  This  scheme  is  based  on  equations  (22)  and  (23), 
but  with  notable  changes  so  as  to  allow  numerical  evaluation  to 
proceed  differently.  The  wing  is  divided  into  a  gridwork 
pattern  as  shown  in  figure  2,  guidelines  for  laying  out  this 
grid  pattern  will  be  brought  out  later.  The  loading  on  the 
wing  is  considered  to  be  given  in  terms  of  concentrated  loads 
rather  than  in  terms  of  distributed  pressure  functions;  this 
is  one  of  the  key  concepts  of  the  approach.  The  concentrated 
loads  are  located  at  the  grid  intersection  points.  Then  by 
means  of  a  "modified"  kernel  function,  downwash  values  are 
derived  at  each  mid-station  point,  illustrated  by  the  point 
labelled  A  in  figure  2;  the  logic  behind  the  choice  of  the  mid- 
station  is  also  brought  out  subsequently.  The  result  is  a 
matrix  equation  which  relates  the  loads  P  in  terms  of  the 
downwash  w  .  This  matrix  equation  is  then  inverted,  giving 
the  desired  end  result  of  P  in  terms  of  w  . 


The  second  key  point  in  the  development  of  this  scheme  is 
the  manner  in  wnich  the  kernel  function  is  handled.  The  strong 
singularities  of  the  kernel  function  along  the  x-axis  aft  of  the 
unit  load  preclude  its  direct  use  in  this  scheme.  An  appreciation 
of  the  type  of  singularities  involved  can  be  obtained  by  picturing 
the  downwash  that  is  associated  with  a  very  narrow  horseshoe  vor¬ 
tex  in  incompressible  flow,  as  shown  in  figure  3(a)  .  The  down- 
wash  pattern  is  singular  in  both  the  x  and  y  directions  at 
the  origin,  and  strongly  so  in  the  y-direction  along  the  horse¬ 
shoe.  The  picture  becomes  all  the  more  dramatic  if  we  envision 
the  situation  as  the  width  of  the  horseshoe  becomes  zero.  To 
remove  the  difficulty  of  dealing  with  these  singularities  we 
modify  the  kernel  function  by  averaging  the  velocities  over 
equally  spaced  intervals  in  the  y-direction,  specifically  we 
write 

,  X 

y9,+  S 

K(x,y0)  =  |  J  K(x,y)dy  (25) 
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where  X  is  the  grid  spacing  in  the  y-direction.  With  this 
simple  modification,  the  downwash  function  changes  to  the  type 
illustrated  in  figure  3(b).  Slice  a  typifies,  for  example, 
the  downwash  in  the  y-direction  ahead  of  the  unit  load;  slice 
b,  a  station  x  aft  of  the  load.  This  averaging  concept  is  based 
on  the  idea  that  in  any  real  situation,  the  wing  merely  senses 
an  average  effect.  Thus,  by  this  averaging  technique,  the 
singularities  in  the  y-direction  are  removed  from  further  con¬ 
sideration.  Note,  a  singularity  remains  in  the  x-direction  at 
the  unit  load,  but  it  does  not  concern  us  because  we  avoid  any 
downwash  consideration  at  this  point  (if  desired,  however,  we 
can  also  consider  an  averaging  out  process  in  the  x-direction 
to  eliminate  this  singularity) . 

In  terms  of  equation  (25),  we  summarize  this  scheme  in 
terms  of  a  matrix  equation  as  follows: 


or 

|w|  =  [K] | p|  (26) 

where  K-^  denotes  simply  the  "average"  value  over  the  interval 
X  in  the  y-direction  of  the  downwash  at  station  m  due  to  a 
load  at  station  n  .  Inversion  of  equation  (26)  gives  the  de¬ 
sired  end  results. 

I  p|  =  [K]  |  w|  (27) 

There  are  several  noteworthy  points  to  bring  out  with  res¬ 
pect  to  the  use  of  equation  (26): 

(1)  No  load  or  pressure  distributions  have  to  be  assumed. 


J 


(2)  It  Is  not  necessary  to  force  any  load  or  pressure  dis¬ 
tribution  to  be  zero  at  the  trailing  edge  (the  condi¬ 
tion  of  zero  pressure  jump  across  the  trailing  wake 
sheet  will  be  satisfied  automatically  because  this 
condition  is  satisfied  for  each  of  the  concentrated 
loads) . 

(3)  The  question  of  how  many  collocation  points  should  be 
used,  and  where  to  locate  them  does  not  enter. 

(4)  Quite  significantly,  this  procedure  affords  a  way  for 
treating  control  surfaces  as  well;  it  is  only  necessary 
to  distribute  concentrated  loads  on  the  control  surfaces 
in  the  same  manner  as  on  the  wing,  and  to  process  these 
simultaneously  with  the  wing  loads. 

(5)  The  averaging  technique  eliminates,  as  mentioned,  any 
difficulties  with  singularities. 

(6)  The  procedure  is  very  routine  and  systematic;  on 
evaluation  and  tabulation,  the  K  values  become 
universal;  che  establishment  of  the  K  matrix  for 
any  case  becomes  simply  a  procedure  of  look  up,  as 
in  the  use  of  trigonometry  tables. 

(7)  The  K  matrix  has  certain  symmetrical  properties 
and  many  elements  are  identical.  For  example, 

=  Kpq  where  p  and  q  have  the  same  position 

relative 4t6  one  another  as  do  m  and  n  ;  also 
Kjnn  =  Kpn  where  m  and  p  denote  stations  of 

equal  distance  to  the  right  and  to  the  left  of  n 
(the  y-direction) . 

(8)  The  scheme  applies  to  any  planform. 

Scheme  II.-  In  this  scheme  the  essential  building  block 
is  the  dipole source  rather  than  the  kernel  function,  and  ex¬ 
plicit  consideration  is  given  to  the  wake .  Steps  are  illus¬ 
trated  in  figure  4  and  are  basically  as  follows.  First,  we 
divide  the  wing  into  a  grid  pattern  as  in  Scheme  I,  and  affix 
a  system  of  concentrated  loads.  Next  we  apply  equation  (14) 
along  each  chordwise  gridline  to  establish  the  strength  of  the 
associated  dipole  line  sources  for  velocity  potential  that  lie 
in  the  chord  plane  of  the  wing.  These  line  sources  are  func¬ 
tionally  continuous  between  the  concentrated  loads,  and  jump  at 
each  load  location.  We  note  here  the  reminder  that  with  respect 
to  equation  (14),  0o(x, y)  =  -  20+(x,y,o)  ,  and  p0(*,y) 

=  -  2p+(x,y,o).  Essentially,  envision  that  we  derive  a  system 
of  dipole  line  sources  which  when  operated  -upon  by  equation  (11) 
yields  back  the  concentrated  load  system  in  the  limit  z  -►  o  . 
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We  then  replace  this  line  source  distribution  by  a  system  of 
"equivalent"  discrete  dipoles,  located  at  the  mid-station 
points;  the  strength  of  each  dipole  is  given  by  €#0(x,y)  f 
where  0O  is  the  strength  of  the  line  source  at  the  mid¬ 
station.  Finally,  we  find  the  total  downwash  at  the  mid¬ 
station  grid  points  on  the  wing  due  to  all  the  discrete  di¬ 
poles  -  those  in  the  wake  as  well  as  in  the  wing  planform. 

A  matrix  equation  similar  to  equation  (26)  results  which 
gives  the  downwash  in  terms  of  the  concentrated  loads;  this 
equation  is  then  inverted,  as  in  Scheme  I.  The  application  of 
this  sequence  of  thought  may  seem  cumbersome,  but  in  practice 
is  quite  simple  and  systematic . 

In  using  equation  (14)  to  establish  the  line  source  dis¬ 
tribution,  we  treat  each  concentrated  load  in  the  nature  of 
a  Dirac  function.  For  a  single  concentrated  unit  load,  located 
a  distance  xE  downstream  from  a  reference  origin,  the  line 
source  distribution  for  potential,  as  established  by  equation 
(14),  is 

+  1  lT"(xr  "  x) 

</>0(x,y)  =  -  20  (x,y,o)  =  -  — e  l(x  -  x  )  (28) 


where  l(x  -  xE)  denotes  a  unit  step  function  at  x  =  x~  . 

The  sequence  or  using  this  equation  for  a  row  of  concentrated 
loads  and  the  means  for  establishing  the  system  of  equivalent 
concentrated  potential  dipoles  is  illustrated  In  figure  5. 


Two  versions  for  treating  the  downwash  from  the  individual 
equivalent  dipoles  bear  investigation.  These  are  depicted  in 
figure  6.  In  the  first  version  we  average  the  downwash  over 
equally  spaced  Intervals  in  the  y-direction  only,  while  In 
the  second  version  an  averaging  is  made  in  both  the  x  and  y 
direction.  This  averaging  removes  the  singularity  problem 
associated  with  the  dipoles.  These  two  averaging  techniques 
are  given  by  the  following  equations  in  terms  of  the  potential 
for  a  dipole  of  unit  strength. 
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Examples  later  on  will  bring  out  the  effect  of  averaging  or 
not  averaging  in  the  x-direction. 

It  should  be  noted  that  the  equation  relating  w  and  P 
as  obtained  by  this  method  should  be  essentially  the  same  as 
equation  (26).  The  only  difference  might  be  that  of  a  slight 
variation  in  the  individual  values  of  the  coefficients  , 

due  to  the  use  of  different  numerical  procedures. 

All  the  points  relative  to  the  utility  of  Scheme  I  apply 
to  this  scheme  also,  but  some  further  points  may  be  made.  The 
fact  that  the  wake  must  be  considered  explicitly  may  at  first 
make  the  scheme  appear  unwieldy.  This  is  really  not  so,  how¬ 
ever;  the  influence  of  downstream  wake  points  decays  very 
rapidly,  and  thus  from  a  practical  point  of  view  we  can  ignore 
the  wake  beyond  a  certain  number  of  chords  distance;  that  is, 
we  simply  terminate  the  wake.  An  example  presented  later  will 
bring  out  the  effect  of  wake  termination.  With  respect  to 
giving  explicit  consideration  to  the  wake  we  note  also  that, 
very  significantly,  we  can  turn  this  fact  to  our  advantage . 
Since  we  have  effectively  established  the  complete  system  of 
potential  dipoles  for  our  wing  system,  we  can  calculate  all 
of  the  induced  velocities  u,v,w  ,  anywhere  in  the  field  as 
well  as  on  the  wing.  This  fact  means  that  we  are  not  res¬ 
tricted  to  planar  lifting  surfaces,  and  that  the  interaction 
effects  of  various  other  lifting  surfaces  in  the  vicinity  of 
the  wing  may  be  handled  as  well.  The  big  additional  advantage 
that  is  offered  by  this  scheme  is  therefore  the  fact  that 
various  configurations  such  as  planar  wings,  wings  with  jumps, 
bi-  or  multi-wing  configuration,  T-tails,  can  be  analyzed. 


Scheme  III.-  This  scheme  represents  a  combination  of 
schemes  1  and  II.  The  key  idea  in  this  scheme  is  the  recogni¬ 
tion  that  the  effect  of  entire  wake  may  be  expressed  in  terms 
of  a  concentrated  load  just  off  the  wing  at  the  trailing  edge . 
This  concentrated  load  is  not  part  of  the  load  on  the  wing,  but 
the  downwash.  it  creates  on  the  wing  is  precisely  that  contributed 
by  the  wake.  It  seems  strange  that  this  fact  has  not  been  re¬ 
cognized  or  utilized  in  previous  analyses.  We  will  deal  only 
with  the  wake  treatment  in  this  section,  since  otherwise  the 
scheme  is  the  same  as  in  scheme  II. 


We  consider  the  situation  that  develops  along  one  of  the 
chord  gridlines  and  especially  the  extension  into  the  wake 
region.  If  we  denote  by  m  the  number  of  concentrated  loads 
that  are  chosen  along  a  chord,  then,  by  the  process  that  is 
illustrated  in  figure  5,  we  find  that  the  intensity  per  unit 
length  of  the  wake  line  source  is  given  by 
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With  reference  to  equation  (28),  this  equation  may  he  written 
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Equation  (31)  is,  however,  precisely  the  line  source  intensity 
of  the  wake  -chat  would  result  from  a  concentrated  load  of 
magnitude  Pe  acting  at  station  xe  .  This  concentrated  load 
acting  just  off  the  trailing  edge  can  therefore  be  used  in  place 
of  the  wake  to  deduce  the  downwash  velocities  that  are  due  to 
the  wake.  For  the  case  shown  in  figure  5,  where  m  =  4  ,  the 
concentrated  load  Pe  would  be  located  midway  between  0O4  and 
0oc  ,  and  would  replace  all  the  wake  to  the  right  of  this 
location.  The  equivalent  wake  load  is  handled  by  use  of  the 
kernel  function  as  in  scheme  I.  It  is  to  be  understood  that  an 
equivalent  concentrated  wake  load  Is  located  at  the  aft  end  of 
each  grid  chord  line,  and  that  of  course  the  superimposed  down- 
wash  effect  of  all  these  loads  must  be  considered. 

In  summary  for  this  scheme,  over  the  wing  region  we  use  the 
equivalent  concentrated  dipole  concept  of  scheme  II,  and  in 
place  of  the  wake  make  use  of  the  equivalent  concentrated  load 
given  by  equation  (32).  It  is  significant  to  note  that  we  deal 
only  with  the  downwash  velocities  ahead  of  this  equivalent  con¬ 
centrated  wake  load  in  treating  the  wing.  There  are  no  singu¬ 
larities  in  this  region,  as  the  kernel  function  so  Indicates 
(except  at  the  concentrated  load,  which  region  does  not  concern 
us)  and  hence  a  singularity  problem  associated  with  the  concen¬ 
trated  wake  forces  is  not  involved. 

Scheme  IV.-  We  present  here  only  a  brief  discussion  of  a 
4th  possible  scheme.  The  basic  idea  of  the  scheme  is  to  start 
with  the  strength  of  the  potential  dipole  sheet  that  is  present 
on  the  wing  and  in  the  wake,  and  to  represent  this  strength 
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bv  convenient  distributions,  say  in  the  form  of  continuous 
functions,  step  functions,  or  concentrated  sources.  The  down 
wash  and  the  loading  would  follow  from  this  distribution  by 
equations  of  the  type 


|w|  =  LA-jjl^l 

|p|  =  LA21 1 0I 


The  elimination  of  from  these  equations  would  in  turn 

establish  a  relation  between  P  and  w  ,  similar  to  those 
derived  in  the  previous  schemes. 


The  details  of  this  scheme  are  left  for  further  develop¬ 
ment  A  derivation  of  the  loading  that  is  associated  with  a 
concentrated  potential  dipole  is  considered  to  be  of  interest 
however,  and  will  be  given  here.  Assume  a  rectangular  strength 
distribution  in  the  x-direction  as  shown  in  the  following 
sketch 


In  the  limit  as  e  -  o  this  distribution  becomes 
dipole  of  unit  strength.  Application  of  equation 
distribution  yields  the  following  force  diagrams 


,a  concentrated 
(11)  to  this 


Thus  we  see  that  a  potential  dipole  of  unit  strength 
upon  by  a  force  of  magnitude  -top  a  °f 

-n  U  These  results,  although  established  here  in 
elementary  way,  may  also  be  derived  by  a  more  formal, 
that  is,  apply  equation  (11)  to  equation  (16),  then 
the  pressures  and  the  moment  of  the  pressures  about 
by  integration  along  the  entire  x-axis,  assuming  z 


is  acted 
magnitude 
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SECTION  IV 
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VERTICAL  VELOCITIES  FOR  UNIT  SOURCE  ELEMENTS 


Potential  dipole.-  The  vertical  velocities  in  the  x  -  y 
plane  for  unit  potential  dipoles,  such  as  those  given  by 
equations(l6)  through  (18),  can  be  established  in  a  very 
simple  unconventional  way  by  a  technique  described  in  this 
section.  We  show  the  method  here,  not  so  much  because  it 
eases  the  treatment  of  dipoles,  but  because  it  illustrates 
the  technique  to  be  used  in  the  next  section  for  treating 
unit  loads,  for  which  cases  the  method  leads  to  a  marked 
simplification  over  conventional  approaches. 

The  formal  way  to  determine  the  downwash  due  to  a 
source  potential  is  to  use  equation  (7c)  and  then  to  evaluate 
the  results  for  z  =  o.  We  can  simplify  the  process  by  expressing 
equation  (7c)  in  difference  form  as  follows 


w 


_  0(z)  -  <p(p) 
z 
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In  all  the  cases  of  concern  to  us  in  this  paper,  however, 
0(o)  =  o  for  x  ^  o  and  y  ^  o;  thus,  for  all  regions  except 
at  the  origin,  w  is  given  simply  by  the  potential  equation 
with  the  z  in  the  numerator  removed,  and  with  z  set  equal 
to  zero  elsewhere  in  the  equation.  At  a  glance,  equations 
(16),  (17)  and  (18)  yield,  for  example: 


For  equation  (l6), 
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For  equation  (17), 
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For  equation  (18), 


These  equations  apply  to  all  points  except  the  origin. 
The  flow  through  the  origin,  where  a  singularity  is  involved, 
is  handled  conveniently  by  use  of  the  finite  grid  concept; 
specifically,  we  obviate  the  necessity  for  considering  the 
singularity  in  detail  by  establishing  the  average  velocity 
through  a  grid  area  which  surrounds  the  origin.  As  an 
example,  the  central  "standpipe"  of  the  velocity  profiles 
shown  in  figure  6  depicts  this  average  velocity.  In  effect 
we  have  replaced  the  point  source  dipole  by  an  equivalent 
dipole  of  finite  size,  very  much  like  the  propeller  dipole 
system  shown  in  the  following  sketch 


The  magnitude  of  the  average  velocity  through  the  central 
grid  area  is  simply  the  sum  of  velocities  through  all  other 
grid  areas,  a  fact  which  follows  from  an  overall  continuity 
condition;  stated  in  a  physical  sense,  we  note  that  all  the 
fluid  that  goes  up  must  come  down. 

For  application  purposes  we  need  to  evaluate  the  velocity 
profiles  for  the  dipoles  of  interest.  Two  will  be  considered 
here.  For  the  case  of  a  subsonic  dipole,  numerical  integration 
of  equation  (J>5)  is  required.  Specifically,  patterned  after 
equation  (30),  we  write  the  average  value  of  velocity  over  a 
grid  area  element  as 
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coordinates  of  the  center  nr  _ n 

we  find  the  average  velocity  over^thn  '  Frora  this  e9uation 
the  origin  to  be  y  over  the  grid  area  surrounding 


"o  -  -Z»n 


n 


(38b) 


flow, ^velocity^rofiles^ayHbe^found  b°r+h"d  lncoraPressible 

of  equation  (36).  The  averse  v!?  by  the  exact  integration 
having  a  a  , 


w  = 
n 


1 

^7T€ 


/ 


x  -+4- 

on  2 


n 


€ 

T 


1 

7 


dx 


7T62  4n2  -  i 


(39a) 


The  center  velocity  w  is 


then 


00 
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(40b) 


00 


A  comparison  between  the  use  of  equations  (39)  and  equations  (40) 
will  be  made  later. 

Unit  load.-  Three  approaches  may  be  advanced  for  determining 
the  downwash  velocities  that  are  produced  by  a  unit  load: 

1. )  Direct  formal  evaluation  by  equation  (22),  as  in 

reference  1. 

2. )  An  approach  which  is  based  on  Scheme  II,  and  which 

makes  use  of  equation  (28);  the  line  source  distribution 
of  potential  dipoles  due  to  the  unit  load  and  given  by 
this  equation  is  replaced  by  "equivalent"  concentrated 
dipoles.  The  velocities  due  to  all  these  concentrated 
potential  dipoles  are  then  superimposed  to  give  the 
downwash  due  to  the  unit  load. 

3. )  An  approach  which  makes  use  of  equation  (13)  and 

equation  (34).  As  in  the  determination  of  equations 
(35) ,  (36),  and  (37),  this  approach  greatly  simplified 
the  determination  of  the  downwash  values  due  to  a  unit 
load,  and  in  the  course  obviates  any  problem  with 
singularities.  The  remainder  of  this  section  is  devoted 
mainly  to  the  development  of  the  basic  notions  of  this 
approach . 

By  means  of  figure  7  we  give  a  quick  and  overall  Insight 
to  the  basic  concepts  of  the  third  approach.  A  detailed  consid¬ 
eration  will  then  follow.  Consider  a  unit  load:  by  the  anti¬ 
symmetric  property  of  p  ,  we  regard  this  load  to  be  split 
equally  between  the  upper  and  lower  half  planes,  sketch  (b). 

With  respect  to  the  associated  pressure  dipole  representation,  we 
take  note  of  the  explicit  appearance,  a  "z"  in  equation  (l6)(und 
in  equations  (17)  and  ( 18 ) ) .  The  potential  due  to  the  unit  loud 
is  given  by  equation  (13).  In  the  region  y  ^  o  ,  the  potential 
is  zero  at  z  =  o  ,  and  hence,  from  equation  (34),  the  vertical 
velocity  is  given  by 

w  =  (41) 

Thus,  near  the  x-y  plane  and  for  y  /  o  ,  <t>  is  a  direct  measure 
of  w  ,  see  sketches  (c),  (d),  and  (e).  The  limiting  form  as 
z  -*•  o  follows  easily.  We  let  the  z  of  equation  (4l)  cancel 
the  z  in  the  dipole  pressure  equation,  and  let  z  =  o 
elsewhere;  the  vertical  velocities  are,  thus,  effectively  given 
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as  the  x-integration  of  equation  (l6)  with  the  z  removed. 

For  y  =  o  and  behind  the  unit  load  (x  >  o),  a  0(0)  exists, 
as  given  by  equation  (14).  The  value  Is  easy  to  establish, 
being  related  to  a  simple  Dirac -type  Integration  through  the 
upper  load  of  1/2,  sketch  (g).  Application  of  equation  (34) 
yields  w  ,  sketch  (h);  effectively,  the  right-hand  branch  of 
0(z)  Is  sheared  downward  an  amount  0‘+(o).  The  magnitude  of  0+(o) 
Is  found  through  the  evaluation  of  a  rather  simple  indeterminant 
form  at  x  =  00  .  The  limiting  value  of  w  as  z  -*•  o  for  this 
slice  Is  also  found  easily,  sketch  (i).  It  Is  to  be  noted  that 
ahead  of  the  unit  load  there  Is  no  0(o);  the  evaluation  of  w 
for  this  region  Is,  therefore,  as  simple  as  for  y  /  o  . 

Further,  the  value  of  0(z)  for  positive  x  can  be  shown  to  be 
given  in  terms  of  0(z)  for  negative  x  .  Thus,  the  determin¬ 
ation  of  w  behind  the  unit  load  can  be  expressed  simply  In 
terms  of  the  w  ahead  of  the  load  and  the  constant  0^(0)  .  The 
development  of  the  velocities  by  the  concepts  shown  in  rigure  7 
is  considered  significant;  thus,  we  can  virtually  form  a  visual 
picture  of  hew  the  velocities  develop. 

With  this  background,  let  us  now  consider  a  unit  pressure 
dipole  in  detail.  We  treat  first  the  regions  y  /  o  /  Substi¬ 
tute  equation  (16)  into  equation  (13),  apply  equation  (34),  and 
let  the  remaining  z's  =  o  .  Since  0(0)  =  o  in  this  region  of 
concern,  the  result  for  w  Is  immediate;  thus, 


where  R  =  v  x  +  p  y  The  notation  K  has  been  introduced 
because  the  expression  is  a  form,  here  greatly  simplified,  of 
the  kernel  function.  With  equation  (18)  and  co  =  o  ,  the  result 
is  simply 


For  use  in  the  grid  system  scr.emes  advanced  in  this  report, 
we  desire  average  values  over  equally  spaced  y-intervals,  in 
accordance  with  equation  (25).  We  have,  therefore,  the  following 
expressions  defining  the  modified  kernel  function  for  regions 
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y  t  O  : 

For  a  subsonic  point  load,  cd  /  o  , 


and  for  a  point  load  with  M  =  o  ,  cu  =  o  , 


K  =  w  = 


For  our  grid  system,  we  use  y0  =  nX  ,  with  n  =  1,2,3>... 
Numerical  evaluation  of  equation  (44)  to  give  tabulated  results 
for  universal  application,  therefore,  is  needed. 

For  the  region  along  and  including  y  =  o  ,  it  is  convenient 
to  proceed  in  terms  of  a  general  notation.  Let  equation  (13)  be 
represented  by 

lOjX 

0(x,y,z)  =  e  l(x,y,z)  (46) 
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where 


px  iu>£/U 

I  =  -  /  e  p(£,y,z)  df  (47) 

J-00 


\ 

An  examination  of  the  Integral  I  in  conjunction  with  the 
pressure  dipole  expressions  given  by  equations  (l6)  through 
(18)  reveals  that  the  integrand  may  be  expressed  in  terms  of 
even  and  odd  functions  of  x  as  follows: 

loog, 

-e  U  p(£,y,z)  =  f(U  +  g(0 


where  f  is  an  even  function  and  g  is  odd.  The  integral  I 
thus  may  be  written 


I(x,y,z)  =  F( 

x,y,z)  +  G(x, 

y,z)  (48) 

where 

F(x,y,z)  =  J 

f U)  at. 

(49a) 

G(x,y,z)  -  j 

•00 

rx 

1  g(4)  ai 

(49b) 

Because 

f  is  even  and  j 

-00 

5  is  odd,  the  follovjing  properties 

exist : 

l(oo)  =  F(oo) 

P(x1)  =  F(») 

> 

-  P(-xx) 

(50) 

G(xx)  =  G(-x. 
F(-oo)  =  G( -00 

L  ^ 

)  =  G(«)  =  0  j 

in  which  the  argument  refers  to  the  x  values.  With  equation 
(34)  in  mind,  we  also  establish  the  function  0(x,y,o)  ,  the 
intensity  of  the  dipole  line  source  on  the  ■+  side  of  the  x-y  plane 
aft  of  the  unit  load.  This  function,  as  obtained  from  equation 
(14),  or  a,s  given  by  -  1/2  the  value  indicated  by  equation  (28) 
with  xr  =  o  ,.  may  be  written 
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(51) 


0(x,  y,o) 


_  lccx 

-  -Xre  U  i6„(o)l(x  -  o) 


PoU 


2  y' 


where  l(x  -  o)  Selection 

erection,  _ “JuMquaMons  ( 3*0 ,  (46) ,  ( 48 )  and  ( 51 )  /the 
vertical  velocity  may  be  written  as 

T?  +  G  -  i  6  (o)l(x  -  o) 

w  ,  HjjiZiZl  -  e  U  - - ^ -  (52) 


-ass "• 

x-axis;  thus 


w(x,z)  = 


Icjox 

1  ‘  u  1 

FU  e  z 


P(x,z)  +  G(x,z )  -  ^  !(x  -  o) 


(53) 


where 


X/2 

F(x,z)  =  x  /  F<x’y'z)  dy 
-X/2 


(54a) 


X/2 


G(x,z)  =  x  /  G<x'y'z)  dy 
-X/2 


(54b) 


We  note  that  the 
x  values,  say  x 


6v(o)  function  has  been  absorbed. 
=  -  x^  ,  equation  (53)  becomes 


For  negative 


Wt-x-^z) 


iCJUX^ 


(55) 


For  positive  x  ,  say  x  -  Xi 
used  to  rewrite  equation  (5d) 


,  equations  (50)  an^  (54)  can  be 
in  the  form 


^(xx,z) 


F(-xx,z)  +  G(-x1,z) 


(56) 
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thus  expressing  w  for  plus  x  In  terms  of  w  for  negative  x  . 

We  now  consider  the  limiting  value  of  these  expressions  as 
z  -*>  o  .  By_the  process  Indicated  earlier  In  this  section,  the 
functions  F/z  and  G/z  reduce  easily  to  the  following 


A(-x)  = 


■rtrfi  0 

-X/2  -oo  R 


103R\  ap 


cos  dy  d£ 

ar 


(57a) 


B(-x)  = 


G(  -x  ,z 


ion 

■  &  r  r  i  o + ^ )e  aps  sin  t* dy  a?  (57b) 


-X/2  -o 


It  is  noted  that  these  equations  are  contained  in  equation  (44) 
as  demonstrated  by  the  following  steps:  suppress  the  factor 

_1  e-lox/U  ^  ]_et  y  =  o  ,  and  split  the  integral  into  a  cosine 
pQu  o 

and  a  sine  part. 

The  remaining  terms  of  equation  (56)  are  combined  to  form 
the  constant 


C  =  I  !(«,)  -  ^ 


(58) 


By  equations  (47)  and  (16),  and  with  a  y-integration,  the  defin¬ 
ition  of  l(°°)  is 


i(«)  -  t 


X/2  »  lffii  2  . 

/  /  e  u  £X(i  ♦  *s§\ 

-X/2  4^r3  V 


-M£  +R 


d£  dy 


/o  2  2  2 

x  +  p  (y  +  z  )  .  The  £ -Integration  can  be  made  by 


means  of  transform  917.8  in  reference 


yielding  the  result 
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X  (-)-xJ  5F  7=F=j  ^  "1  VU  v  J 
-x/2  v  y  +  z 

■i  00  /v2  +  as  will  be  the  case  in  practice , 

For  small  values  of  y  v  y  +  z  >  as  Wilx 

this  expression  reduces  to 

X/2 


I(°°) 


I  F  J£—  _ J  _ _ _ 

"  x  J  "2  “ 

-X/2  1 


— .  .  1  .  -1  X 

I(co)  =  £  tan  22- 


With  this  equation,  C  becomes 


c  “  ikk  (tan  1  ^  ' 


or  in  the  limit  as  z  -►  o 


C  =  - 


(6la) 


(6lb) 


_  find  the  following  equations  apply  for  the  average 
In  summary  }  we  find  i j  4-v, o  +■  f>entpr,ftd  over  y  —  o  i 


velocity  along  the  grid  band  that  la  centered  over  y  -  o 

iU2Cl 

K  =  5(-xro)  =  A  e_!r  [AA>  *  ’  x  =  -X1 


(62a) 


ia«1 

I  -  S(*i-o)  -  p^^[-A(-xl)  +  B(-xl)'^J’  X  =  Xl  (62b 

modifleT^ernel  functlonho^the^presen/repor/and^re^for  use 
in  equation  (26)  of  Scheme  I. 

It  is  noted  that  equations  analogous  to  equations  (62) 
apply  also  for  the  case  of  2-d  flow;  they  are  somewhat  simpler 


and  differ  In  the  following  three  specific  ways:  (1)  the 
y-integration  leading  to  equation  (53)  is,  of  course,  not 
involved;  (2)  the  factor  1/2  would  appear  in  place  of  the 
factor  1/2X  in  equation  (53);  and  (3)  the  constant  C  ,  as 
given  by  equation  (6la)  would  be  different.  A  later  example 
will  illustrate  these  points. 


SECTION  V 


EXAMPLE  SPECIFIC  CASES  FOR  K 


Case  1;  2-dimensional  incompressible  steady  flow  (cd  =  o).- 
To  bring  out  some  interesting  and  instructive  facts,  determination 
of  the  vertical  velocity  produced  by  a  unit  load  will  be  illus¬ 
trated  by  five  different  techniques  mentioned  in  this  paper. 

Solution  1: 


Use  is  made  here  of  equations  analogous  to  equations  (57)  and 
(62) .  For  this  case  we  find 

A<x)  -  -  2^ 

B(x)  =  0 
C  =  0 


From  equations  (62),  these  lead  to  a  single  expression  for  w 
given  by 


w 


-1 

£ttP0Ux 


(63) 


which  agrees  with  the  known  exact  solution  to  the  problem. 
Solution  2; 


This  solution  is  based  on  Scheme  II,  and  is  presented  with  the 
aid  of  figure  8.  For  the  unit  load  the  wake  potential  is  given 
simply  by  (see  equation  (28)) 

0  - - i—  1  (x  -  o) 

^o  p  U  '  ' 

*o 

The  strength  of  the  equivalent  concentrated  potential  dipoles 
spaced  at  equal  intervals  €  is  thus 


The  w’s  from  each  of  these  dipoles  are  given  by  equations  (39); 
the  values  are  tabulated  on  respective  lines  below  the  dipoles. 


The  superimposed  total  w  is  shown  near  the  bottom  of  figure  8; 
it  is  noted  that  the  results  are  the  same  as  given  by  equation 
(63)  j  thus  showing  that  this  method  also  leads  to  the  exact 
answer . 

Solution  3: 


This  solution  proceeds  as  in  the  second  solution,  but  makes  use 
of  equations  (40)  instead  of  equations  (39).  The  results  ob¬ 
tained  by  this  method  are  shown  in  figure  9>  curve  (b),  along 
with  the  exact  results  just  derived,  curve  (a).  We  see  from 
this  comparison  the  effect  of  not  using  average  va'’  ues  of  velo¬ 
cities  over  each  interval  in  the  representation  of  w  for  each 
concentrated  dipole . 

Solution  4: 

In  this  solution  we  determine  w  by  equation  (34),  and  keep  z 
finite  to  show  the  influence  of  not  passing  to  the  limit  z  -►  o  . 
Equations  analogous  to  equations  (55)  and  (56)  are  involved. 

With  equations  (13) ,  (14)  and  (17)  the  result  can  be  shown  to  be 


0(x,z) 


zd£ 

£ +  z 


2 


-  (tarrl  f  + 1) 

1(x  "  °) 

»  -  ^  k  [tan-1  f  +  f  -  *(«  -  °)] 

=  -  f  [tan-1  f  +  f  -  irl(x  -  o)]we  (64) 


where  we  is  the  exact  result  given  by  equation  (63).  The  re¬ 
sults  axe  shown  in  figure  10  for  positive  x  only,  since  the 
results  for  negative  x  would  be  the  same.  We  note  the  in¬ 
teresting  fact  that  as  long  as  £  >  3  >  the  approximate  solution 
(z  finite)  is  essentially  the  saffle  as  the  exact  solution. 
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Solution  5: 

In  this  solution  we  show  the  effect  of  wake  termination.  The 
wake  strength  is  given  by 


*o<x>  =  -  TV  1(x  •  °> 
rO 


With  equation  (17)  we  find  the  field  potential  for  a  finite 
length  wake  to  be 


0(x, z) 


zd£ 


(x  -  I)  +  z‘ 


where  b  represents  the  point  downstream  of  wake  cutoff.  This 
equation  yields 

*(x>z>  -  2^  (tan_1 1  - tan_1  S-6) 


The  vertical  velocity  as  determined  from  this  equation  is 


w  “  llm  If  =  ‘  (}  +  F^bc) 

z  -*  o 


This  equation  shows  that  w  will  be  essentially  the  exact  value 
if  ^  is  kept  small.  A  useful  guide  is  afforded  by  this 

equation j  thus,  the  term  --JS —  shows  that  if  x  is  kept  to 
within  .05  b,  for  example, Dtfieft  errors  in  w  due  to  wake 
termination  will  be  less  than  5  percent.  This  observation,  of 
course,  applies  to  the  2-dimensional  steady  flow  problem.  For 
the  oscillatory  case,  the  errors  in  w  due  to  wake  termination 
should  be  tven  less. 

Case  2:  2-dlmensional  incompressible  flow,  co  7*  o  .-  For 
this  case  equations  corresponding  to  equations  (57)  are 


A(-x) 


1 

27 T 


(66a) 


* 
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d£  =  -  i 


(66b) 


Because  2-dimensional  flow  is  involved,  the  constant  C  for 
this  case  is  different  than  that  given  by  equations  (58)  and  (6l). 
Here,  in  place  of  the  I  in  equation  (56),  we  have 


CO  i si 

r  zeU 

J  (S  +  z2 


1 

2 


-z 


e 


CD 


U 


and  in  place  of  the  factor  we  have  simply  1/2;  the  constant 

C  is  thus  given  by 


1  CD 

“  2U 


Using  equations  (66)  and  this  value  of  C  in  place  of  the - ^ 

term  in  equations  (62),  we  find  rrX 


w(x) 


1 

2ttP0u 


+  1?  +  i  Si(sg£)]}  (67) 


which  agrees  with  the  known  solution  of  this  problem  as  de¬ 
termined  by  other  exact  means;  it  is  interesting  to  note  that 
this  equation  applies  for  both  negative  and  positive  x  . 

Case  3s  3-dimenslonal  incompressible  flow,  cd  =  o  ,-  The 
solution  desired  in  this  case  is  for  the  average  w  ove r 
equally  spaced  y- intervals,  as  would  be  used  in  equation  (26). 
The  equation  for  w  for  this  case  applying  in  the  region 
y  /  o  has  already  been  derived,  equation  (45).  For  y  =  o  , 
we  find 


34 


2n  +  1)*  +  4 

X* 


This  equation  may  be  investigated  quite  readily  for  x  very 
small  and  x  very  large;  the  known  infinite  summation 


w 

1=1 


4n  -  1 


provides  an  aid.  The  results  are  found  to  be  simply 


2tjp  Ux 


,  x  small 


w  =  o 


X  =  00 


These  results  indicate  a  rather  remarkable  fact;  thus,  in 
spite  of  the  fact  that  we  have  represented  the  uniform  line 
source  by  an  infinite  array  of  concentrated  loads,  we  find 
the  solution  to  be  the  exact  solution  for  small  and  large 
x  .  The  use  of  average  values  of  w  over  each  Interval 
evidently  retains  all  the  essential  intelligence  of  the  flow 
field.  Equation  (69)  has  not  been  investigated  here  for  in¬ 
termediate  x  values,  but  it  is  expected  to  yield  exact  re¬ 
sults  in  this  range  also. 
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SECTION  VI 


EXAMPLE  TREATMENT  Off  'wrung 


schemes  oftthiB°pap£r  are  show^l^thi  th®  numerical  grid 
and  because  exact  solutionsare  known* s  section;  for  brevity, 
in  terms  of  a  2-dlmenslonal  alrfSfTSi  ^ady^low?  8lVen 

Example  1; 


grid  spacingPle  decora*  ^t^cho^*6  concen^rated  load.  The 
at  the  quarter-chord  mint-  chord,  c  *  Pl&ce  the  load 
in  figure  S  (a  )°  (our  tasLhf  \  at  3A>  an  shown 

prior  knowledge  of  airfn-n  hov,  these  locations  is  based  on 
tlon  (63)  wl  flnd  that  behaVl0r  oh^“terlstles).  By  equa- 


w 


-L  _ 

Trp  Uc 
ro 


leading  to  the  inverted  form 


L  =  -  irpoUcw 

For  a  constant  angle  of  att-nnir  n  , 

thus  B  attack  a  ,  we  have  w  =  -uo  ,  and 


L  = 


rp/co 


moment  will  Sso^^coJre^t^e^ause  thl  io°dlem‘  The  pltchjLn6 
quarter-chord.  Thus  we  see  the  rather.6  load  was  Placed  at  the 
correct  lift  and  moment  were^btataed usin^o*?16  that  the 
and  a  single  downwash  point.  We  have  Single  load 

absurd  case  because  it  provides  akev  d  tJis  feemlngly 

the  next  examples.  viaes  a  key  guide  on  how  to  handle 

Example  2: 

£  ngurf^.^wuh  e*L°p?eS  Ad°TTh  POlnt*  as  sho““ 

these  loads  and  w  points  in  %fUi?e'  we  have  located 

note  the  locations  with  care  ^iCSar  }£cations  shown; 

pma^a  equations  lo^Sls^se  .«"fS‘W3)  the  appr°- 


W1 

2 

-i  i 

w_ 

= 

i 

2 

-  T  -i 

L2 

Inverted,  these  equations  appear 


L1 

_  ^PqUc 

-3 

-3 

W1 

L2 

8 

1 

-3 

W2 

For  a  constant  angle  of  attack,  for  which  w,  =  w  =  -Ua  ,  we 
find  1  d 

L1  =  \  ^P0U2ca 
l2  =  \  TpQU2ca 

These  values  yield  a  total  lift  and  a  moment  arm  referred  to  the 
leading  edge  as  follows 

L  =  Trpou2ca 

iB  +  5tr_c 

TT-'* 

These  values  also  correspond  to  the  exact  solution. 

Example  3i 

Thia  example  considers  3  loads  and  3  downwash  points,  located 
as  shown  in  figure  12(c).  The  guideline  being  evolved,  and 
which  appears  to  be  remarkably  successful,  is  as  follows: 
divide  the  chord  into  as  many  equal  intervals  as  there  are 
lotdsj  locate  a  load  at  each  1/4  point  of  an  interval,  a  down- 
wash  point  at  each  3/4  point.  The  equations  that  result  for 
this  case  are: 


1 


T 

-15  -10  -15 

W1 

TTp  UC 

L2 

64 

6  -12  -10 

w2 

L3 

1  6  -13 

w3 

For  a  constant  angle  of  attack,  with  w-,  =  wp  =  w_  =  -Ua  ,  we 
find  3 

I.,  =  ^  TrpoU2ca 

*  |  1rPou2°Ct 

l3  ”  H  tPou2c0 

which  yield  a  total  lift  and  a  moment  arm  relative  to  the  leading 
edge  of 

L  =  irpoU2ca 
c 

€  =  S 

Again  the  exact  results  are  obtained.  By  induction,  we  may  expect 
an  exact  solution  for  all  Bmaller  intervals.  If  the  loads  and 
downwash  points  are  located  by  the  guideline  stated  at  the  begin¬ 
ning  of  this  example.  From  these  three  examples  we  wish  also  to 
make  the  following  observation  thought  significant.  In  all  the 
three  cases,  downwash  points  have  been  chosen  midway  between  the 
concentrated  loads  and  in  an  aft  direction,  with  surprisingly 
successful  end  results.  No  consideration  has  had  to  be  made  about 
flow  or  loading  in  the  leading  edge  region.  Essentially,  we  are 
saying  that  with  this  technique,  we  care  not  a  wit  about  the  math¬ 
ematical  character  of  the  loading  at  the  leading  edge,  nor  whether 
w  is  satisfied  in  this  region.  Apparently  flow  is  controlled  pre¬ 
dominantly  by  the  aft  regions  of  the  airfoil.  Perhaps  this  is  a 
reflection  of  the  Kutta  condition  at  the  trailing  edge,  which 
governs  the  strength  of  the  entire  flow. 

Example  4: 

For  this  example.  Scheme  II  was  applied  to  the  airfoil  case 
treated  in  Example  3.  A  detailed  listing  of  the  solution  is 
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not  considered  necessary,  but  solution  proceeds  in  the  manner 
depicted  by  figures  5  and  8.  Again  considering  the  constant 
angle  of  attack  case,  it  was  found  that  exact  results  for  lift 
and  moment  were  also  obtained  by  this  approach.  A  main  point 
of  the  example  is  the  reminder  that,  since  the  complete  poten¬ 
tial  system  is  constructed  by  this  approach,  it  is  possible 
to  establish  all  components  of  flow  velocity  anywhere  in  the 
field,  such  as  in  the  vicinity  of  a  tail. 
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SECTION  VII 


REMARKS  ON  GRID  LAYOUT 


The  concepts  advanced  in  this  paper  have  not  been 
exercised  sufficiently  to  permit  definite  guidelines  to  be 
stated  on  wing  grid  layout.  Some  suggestions  on  how  to  lay 
out  the  grid  system  can  be  made,  however.  It  is  felt  that 
a  matrix  somewhere  between  a  4x4  matrix  (four  chordwise 
intervals  by  four  semispan  intervals )  and  a  10  x  10  matrix 
should  adequately  handle  most  lifting  surfaces.  The  thought 
that  the  matrix  should  be  square  is  not  to  be  inferred  from 
these  numbers.  Suggested  layouts  are  shown  in  figure  13. 

For  untapered  straight  or  swept  wings  and  for  deltas  with  a 
straight  trailing  edge,  a  simple  rectangular  grid  pattern 
can  be  used,  figure  13(a)  and  13(b)*  For  deltas  with  a 
swept  forward  trailing  edge,  an  array  consisting  of  rec¬ 
tangular  blocks  of  one  size  in  the  leading  edge  region  and 
rectangular  blocks  of  another  size  in  the  trailing  edge 
region  might  be  employed,  figure  13(b);  these  smaller  trail¬ 
ing  edge  blocks  represent  the  type  that  might  also  be  used 
in  the  treatment  of  control  surfaces.  An  interesting  ques¬ 
tion  arises  at  this  point.  Can  intervals  of  different 
lengths  be  mixed?  As  a  means  for  shedding  some  light  on 
this  question,  the  wing  system  shown  in  the  following  sketch 
was  analyzed. 


The  airfoil  was  arbitrarily  divided  into  two  intervals  -  the 
front  interval  being  4/5  c  long,  the  aft  1/5  c  long.  A  load 
and  a  downwash  point  were  located  in  each  interval  according 
to  the  rule  of  thumb  used  in  the  example  airfoil  cases;  spec¬ 
ifically,  the  load  is  placed  at  the  1/4  point  of  the  interval, 
the  downwash  at  the  3/4  point.  The  system  was  analyzed,  as 
in  the  example  cases,  with  the  following  results: 

Lx  =  .923  ttp0U2ccx 

L2  =  .077  irp0U2ca 

L  =  irp0U2ca 

e  =  c/4 
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Amazingly,  the  exact  answer  is  found  again.  (There  appears 
to  be  some  magic  associated  with  the  1/4  and  3/4  points.) 
Evidently,  then,  mixed  interval  size  can  be  used  to  good 
advantage . 

Figure  13(c)  shows  thoughts  on  tapered  wings.  We  note 
that  the  blocks  are  of  varying  aspect  ratios  in  this  constant 
y-interval  treatment.  A  variation  of  this  scheme  is  shown  in 
figure  13(d)  wherein  the  y-interval  is  also  varied  to  yield 
blocks  of  constant  aspect  ratio.  The  possible  disadvantage 
of  the  layouts  of  figures  13(c)  and  13(d)  is  that  unequal 
intervals  and  a  sheared  grid  system  are  involved,  thereby 
increasing  the  task  of  evaluating  the  kernel  function;  this 
point,  of  course,  does  not  become  of  too  great  a  concern  if 
the  kernel  function  values  are  tabulated  at  frequent  x  and 
y  intervals.  In  spite  of  the  fact  that  the  load  and  down- 
wash  points  are  not  located  in  a  rectangular  array  fashion, 
the  layout  scheme  of  figure  13  (d )  may  turn  out  to  be  the  pre¬ 
ferred  scheme.  Of  all  the  schemes,  it  is  the  least  arbitrary, 
and  it  applies  to  all  wings  regardless  of  planform.  The 
scheme  also  tends  to  distribute  and  group  the  loads  more 
towards  the  tip  sections,  which  probably  are  the  most  important 
in  controlling  aeroelastic  behavior. 

Figure  13(e)  shows  a  possible  gridwork  for  a  trapezoidal- 
type  wing;  here,  a  fairly  fine  rectangular  gridwork  Is  used. 

Two  additional  points  may  be  brought  out  In  discussing  this 
arrangement.  One  pertains  to  gridwork  size  In  the  Inboard 
region.  If  small  blocks  are  used  throughout,  then  many  chord- 
wise  blocks  are  found  In  this  region.  It,  thus,  may  be  worth¬ 
while  to  Increase  the  size  of  these  inboard  blocks  as  a  means 
for  decreasing  the  total  number.  In  the  example  shown,  58 
blocks  would  be  involved  if  the  same  small  size  were  used 
throughout.  By  only  doubling  the  length  of  the  inboard  blocks, 
the  number  has  been  reduced  to  38.  This  Idea  should  especially 
apply  if  the  wing  is  a  cantilever,  since  the  Inboard  region 
probably  has  little  influence  on  determining  the  flow,  and 
therefore  should  not  require  as  detailed  a  treatment  as  for 
the  tip  region.  A  second  point  is  that  with  the  use  of  small 
blocks  it  probably  isn't  necessary  to  average  the  velocity  in 
the  y-dlrection  for  blocks  remote  to  a  load  (as  in  equation 
(44)).  The  use  of  the  w  at  the  center  of  such  blocks  is 
probably  accurate  enough,  thus  lessening  the  chore  of  deter¬ 
mining  the  kernel  function.  The  effect  of  averaging  or  not 
averaging  can  be  Illustrated  by  equation  (39a);  the  -1  in 
the  denominator  is  a  result  of  the  averaging  process,  and 
without  the  -1  the  equation  would  yield  the  velocity  at  the 
center  of  the  interval.  It  can  be  seen  that  for  n  as  small 
as  4  ,  referring  to  a  point  just  4  intervals  away,  the  neglect 

of  the  -1  changes  the  w-value  by  not  even  2  percent. 


42 


7 


How  large  should  the  Intervals  e  and  X  be?  The 
steady  state  example  cases  Indicate  that  the  interval  may  be 
quite  large  (even  equal  to  the  chord  in  one  case).  For  oscil¬ 
latory  flow,  however,  definite  wavelength  patterns  exist  in 
the  flow  field,  and  these  can  be  expected  to  have  a  bearing 
on  interval  size.  We  consider  these  wavelengths  in  the  following 
as  a  means  for  establishing  an  upper  limit  to  interval  size. 

An  examination  of  equation  (16)  indicates  that  three  frequency 
components  in  the  x-direction  must  be  considered  if  evaluation 
is  made  by  numerical  techniques;  namely. 
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03(1  -  M) 
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TT(T  +  M) 


„  03(1  +  M) 
op2 
oM 

=  TT(T"-'M) 


for  x  positive,  R  =  x 


for  x  negative,  R  =  fx 


The  smallest  wavelengths  are  associated  with  ft,  and  ft. 
and  are  given  by 


X 
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2ir  _  2rrU 

1  17^  03 

2v  2irU(l  -  M) 

3  =  rq  =  — 


To  maintain  a  reasonably  accurate  numerical  procedure,  we 
reason  that  the  interval  e  should  be  something  in  the  order 
of  1/12  or  less  of  the  shortest  wavelength;  (effectively, 
we  are  saying  that  12  equally  spaced  values  or  more  should 
define  one  cycle  of  a  sine  wave  reasonably  well).  Thus 


p  .  2irU 
6  ^  1255 


or 


c  <  • 

whichever  is  the  shorter, 
reduced  frequency,  we  have 


2ttU(1  -  M) 

"120*4 - 

Written  in  terms  of  the  chord  and 
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1  ^  JL. 

c  —  T2lc 


e  x  JL-  1 

c  -  TSE  M 


The  firat  of  theee  equation^governScfo^  t^e 

Expressions  ail  given’ In  the  following  tabulation. 
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.05 

.1 

.2 

.3 

.4 

•  5 


M  =  o  -  .  5 
e/c  =  5-3 
2.62 
1.31 
.882 
.654 
.523 


A 

3-48 

1.74 

.882 

.582 

.436 

.348 


A 

1.31 

.654 

.327 

.196 

.164 

.131 


A  similar  treatment  for  the  y-dlrection  Indicates  only 
a  slngle^frequency  and  the  following  results. 
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These  results  indicate  that  the  wavelei»gth  components  of 
the  flow  are  not  of  prime  concern  in  establishing  interval  size, 
except  in  the  range  of  high  k  and  high  "subsonic  Mach  number 
values.  For  the  lower  k  and  M  values,  other  considerations, 
such  as  the  geometry  of  motion  of  the  airfoil,  would  govern. 

With  respect  to  a  plunging,  pitching,  and  possibly  deforming 
airfoil,  for  example,  the  following  question  is  of  concern. 

How  many  downward  points  should  be  used  to  represent  the  motion? 
From  a  curve-fitting  point  of  view,  it  is  reasoned  that  a  min¬ 
imum  of  4  to  5  points  should  be  used  to  give  adequate  represen¬ 
tation  of  motion  across  the  chord.  Thus  a  minimum  of  4  chord- 
wise  intervals  is  suggested;  by  similar  consideration,  at  least 
4  spanwise  stations  should  be  used.  Taking  into  account  these 
various  thoughts  on  interval  size,  and  until  more  experience  on 
the  practical  application  of  the  grid  scheme  is  available,  we 
might  state  the  following  rule  of  thumb:  for  low  k  and  low  M 
values,  use  5  or  6  chordwise  and  5  or  6  spanwise  Intervals,  and 
for  the  high  k  and  M  ranges,  use  intervals  as  indicated  by 
the  e/c  and  X/c  tabulations  given  earlier  in  this  section. 
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SECTION  VIII 


FREQUENCY  RESPONSE  AND  FLUTTER  DETERMINATION 


One  of  the  features  of  the  schemes  of  this  paper  is  that 
the  equation  relating  P  to  w  is  ideally  suited  for  both 
of  the  basic  types  of  aeroelastic  formulations  that  are 
normally  employed,  specifically,  for  either  the  approach  which 
is  in  terms  of  discrete  masses  and  influence  coefficients  or 
the  approach  which  uses  modal  functions  and  a  Lagrangian  treat¬ 
ment.  By  way  of  example,  we  show  in  the  following  how  the 
frequency  response  functions  of  an  airplane  due  to  sinusoidal 
gust  encounter  may  be  derived;  the  equations  for  flutter  are 
in  turn  automatically  derived. 

In  reading  this  section  the  reader  should  be  aware  that, 
to  simplify  the  writing,  considerable  liberty  has  been  taken 
with  the  matrix  notation;  it  is  abreviated  or  inconsistent  in 
places,  and  matrix  elements  are  presented  in  a  general  sense 
only. 


By  discrete  masses  and  influence  coefficients.-  Let  the 
dynamical  equations  for  sinusoidal  motion  of  the  complete  air¬ 
plane  system  be  represented  by  the  matrix  equation. 

[D]|zs|  =  <o2|m||zs|  +  |PS|  (70) 


where  Z§  is  the  column  matrix  denoting  the  deflection  at  each 
of  the  discrete  mass  points  (which  may  include  the  fuselage  and 
tail  as  well  as  the  wing),  D  is  the  operator  which  leads  to 
total  loading,  m  a  diagonal  mass  matrix,  and  Pq  the  applied 
loading  matrix.  In  the  form  presented,  the  equation  applies 
either  to  a  free-free  system  or  to  a  restrained  system  such  as 
a  cantilever.  For  a  sinusoidal  gust,  equation  (24b)  may  be 
written 

^cux 

w  =  U  |§  +  iouZ  -  e  U  (71) 


where  x  is  the  location  of  the  downwash  point  relative  to 
some  convenient  reference  point  such  as  the  leading  edge  of 
the  root  chord  of  the  wing.  In  general,  the  load  and  de¬ 
flection  points  used  for  the  dynamical  treatment  of  the  system 
may  not  be  the  same  is  those  used  for  establishing  the  aero¬ 
dynamic  loads.  By  suitable  interpolation  formulae,  however, 
these  two  load  and  deflection  systems  may  be  easily  related, 
and  we  suppose  that  the  transformations  from  one  system  to  the 
other  are  given  by 
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A 


|P|  -  (Til  |PSI 

(72a) 

|z8l  =  (t21  lzl 

(72b) 

where 

Combining 

and  To  are  in  general  rectangular  matrices, 
equations  (26),  (70),  (71)  and  (72),  we  find 

_i  — 

K  TX(D  -  o)2m)T2Z  =  (u  ^  +  ioa)  Z  -  e  U 

(73) 

where,  for  simplicity  in  writing,  the  matrix  notation  has  been 
dropped.  Equation  (73)  in  turn  may  be  written 


K  TX(D  -  o&n)^  -  (u  ^  +  i^)  Z  -  -  e 


or 


_i 


DXZ  = 


-  e 


(UK 

U 


(74) 


Inverted,  we  find  the  desired  result  for  the  frequency  response 
function  Z  ,  thus 


i 

Z  =  -  D1'*1  e  U  (75) 

Prom  Z  we  may  determine  the  frequency  response  function  for  I 
other  quantities  of  interest,  such  as  acceleration  or  load  at 
given  points. 

It  is  to  be  noted  that  this  development  was  made  in  a  way 
so  as  to  make  the  inversion  of  the  R  matrix  unnecessary,  see 
equation  (74) j  thus  the  developments  in  this  report  appear  all 
the  more  attractive.  The  Dt  matrix  of  equation  (74)  is  ob¬ 
tained  by  matrix  multiplication  and  addition  of  various  matrices, 
each  one  of  which  is  established  in  a  straightforward  way. 
Solution  for  the  frequency  response  is  thus  reduced  to  a  single 
inversion  of  the  end  result  matrix  Dt  ,  as  shown  by  equation 
(75). 
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The  equation  describing  flutter  of  the  system  follows 
directly  from  equation  (74),  being  simply  this  equation  with 
the  right  hand  side  set  equal  to  zero. 

By  a  modal  function  approach.-  We  assume  that  the  response 
due  to  a  sinusoidal  gust  input  is  expressed  in  terms  of  the 
natural  modes  Z  of  the  airplane  system  according  to  the 
equation  n 


Z  =  a1Z1  +  agZ2  +  a^Z^  +  ... 


(76) 


where  the  a~' s  are  generalized  coordinates.  It  is  well 
established  that  a  Lagrangian  formulation  in  terms  of  equation 
(76)  leads  to  the  following  basic  response  equation  for  an 

Vn  +  VVn  =/pVS  <77) 


where  Mn  is  the  generalized  mass  and  is  the  natural 

frequency  of  the  nth  mode,  and  where  p  is  the  applied 
loading  over  the  surface  S  .  If  equation  (77)  is  applied  to 
each  of  the  modes  considered,  and  if  the  loading  p  is  ex¬ 
pressed  in  terms  of  equivalent  concentrated  loads  at  various 
grid  points,  then,  for  the  sinusoidal  case  of  ,  we  can 
write  the  following  matrix  equation 


K%2  -  “2k]W  ■  M’l  (78) 

where  [X2  -  “2X]  is  a  simple  diagonal  matrix,  aR  a 

column  matrix,  Zn  is  a  rectangular  matrix  built  up  from  row 
matrices  which  express  the  deflection  at  each  grid  point  for 
each  mode  Zn  ,  and  P  is  a  column  matrix  cf  the  applied  con¬ 
centrated  loads  at  each  grid  point.  The  simple  matrix  form  of 
the  generalized  forces  on  the  right  hand  side  of  equation  (78) 
follows  from  equation  (77)  because  of  the  Dirac  function  nature 
of  each  concentrated  load.  We  now  show  two  ways  for  proceeding 
with  the  solution.  One  method  involves  the  inversion  of  the  K 
matrix  and  leads  to  a  solution  in  terms  of  the  values.  The 

other  rep’cesents  a  very  interesting  different  version,  which 
avoids  having  to  invert  the  K  matrix  and  which  leads  to  a 
solution  directly  in  terms  of  the  P  values.  For  both  versions 
we  make  use  of  the  matrix  representation  of  equation  (76)  to  ex¬ 
press  the  deflection  at  each  of  the  chosen  grid  points,  namely 
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z|  -  MaJ 


where  fZ~J  is  a  rectangular  matrix  is  which  each  column  re¬ 
presents  the  deflection  at  the  various  grid  points  in  a  given 
mode . 

In  the  first  version  we  combine  equations  (27)  and  (71), 
and  (79)  to  obtain  the  loading  as 


(79) 


P|  =  [K)"1[u^  +  ia3]fZ^j|an|  -  {R]-1! 


,  OK 

u  I 


(80) 


The  substitution  of  this  equation  into  equation  (78)  leads  to 
the  result 


.(V  -  “2X]  -  LVfn'X  ^  +  i4Zn 


an 


-  -  Unl  [S]_1|e  *  U  | 


“I  -fr- 


or 


1  d 


jL 

[E]|aJ  =  -  LZnltK]-1|e  U  |  (81) 

Solution  for  an  yields  the  frequency  response  functions  for 
ajj  .  With  the  a^s  established,  various  other  frequency 

response  functions,  such  as  acceleration  or  load  at  a  given  • 

point,  follow  readily. 


In  the  second  version  we  proceed  by  solving  equation  (78) 
for  laj  ,  or 


where  solution  follows  simply  because  the  square  matrix  on 
the  left  hand  side  of  equation  (78)  is  diagonal.  The  combina¬ 
tion  of  this  equation  and  equation  (79)  yields 


With  this  equation  and  equation  (71),  we  find 
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I  =  [u  Jr  +  1<“]rzr.,[ 


<%2  ' 


]lZnll 


_■» 


COX 


P  -  e 


1  TT, 


(82) 


We  now  combine  this  equation  with  equation  (26)  by  eliminating 
w  ,  and  find  the  following  equation  in  terms  of  the  P  values 
only 


[[K]  -  [u  Jr  +  i4znj[(“2  -;  ^)M  ]LZnl] 

v  n  '  n 


-1 


Pi  =  -  I  e 


OK 
U  I 


(83) 


or 


[H  3 1  P  |  = 


-  e 


-i 


cux 
U  1 


Inverted,  we  find  the  frequency  response  values  for  P 
directly  as 


P 


(84) 


It  is  to  be  noted  that  the  matrices  leading  to  H  are  all 
very  easy  to  establish,  and  thus  this  version  appears 
especially  attractive. 

For  flutter  considerations  by  this  modal  approach  we 
consider  the  homogeneous  part  of  either  equation  (8l)  or  (83); 
that  is,  we  consider  these  equations  with  the  right  hand  sides 
set  equal  to  zero. 
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SECTION  IX 


CONCLUDING  REMARKS 


Most  of  the  sections  in  this  paper  contain  some  discussion 
of  the  possible  merits  or  advantages  of  the  concepts  developed. 

We  therefore  conclude  with  only  a  few  additional  remarks. 

In  this  paper,  we  have  advanced  the  concept  of  using 
concentrated  loads  to  represent  the  aerodynamic  loading  on  a 
lifting  surface.  Many  variations  of  this  idea  may  be  conceived. 

For  example,  the  distributed  loading  may  be  expressed  in  terms 
of  various  approximations,  such  as  straight  line  or  parabolic 
segments,  with  unknown  values  at  various  grid  points.  Equivalent 
concentrated  loads  can,  in  turn,  be  deduced  from  this  represent¬ 
ation.  These  concentrated  loads  can  then  be  handled  by  the 
procedures  developed  in  the  paper,  with  the  result  that  the 
downwash  values  are  thus  given  in  terms  of  the  distributed  load 
values.  Whether  this  "apparent"  refinement  would  represent  an 
Improvement  or  not,  however,  is  not  known.  Another  idea  that 
may  be  of  good  practical  use  in  application  is  the  following. 
Consider  a  slice  cut  out  of  the  wing  in  the  x-direction  and  of 
width  X  in  the  y-dlrectlon.  Assume  a  pressure  distribution  on 
this  slice  that  is  uniform  in  the  y-dlrection  and  represented 
by  various  pressure  mode  functions  in  the  x-dlrection.  For  each 
of  the  pressure  mode  assumptions,  establish  the  downwash  for  a 
number  of  points  along  the  center  line  of  the  slice  and  along 
lines  in  the  x-dlrectlon  which  are  spaced  at  equal  y-lntervals 
X  ;  these  downwash  values  are  then  used  for  universal  application 
thereafter.  The  application  to  any  wing,  for  example,  would  be 
to  consider  a  series  of  such  slices  placed  side  by  3ide  in  such 
a  way  as  to  cover  the  planform;  essentially,  by  this  representa¬ 
tion,  we  envision  that  the  pressure  distribution  over  the  wing 
is  given  by  smooth  functions  in  the  x-dlrectlon  and  by  a 
succession  of  steps  in  the  y-dlrection.  Solution  by  this 
approach  becomes  simply  a  matter  of  determining  by  an  equation 
analogous  to  equation  (27)  the  magnitude  of  the  x-pressure  modes 
at  each  y  station  so  as  to  yield  the  prescribed  or  desired 
downwash  distribution. 

With  respect  to  limits  on  matrix  size  in  the  practical 
application  of  the  schemes  of  this  paper,  we  may  find  that  there 
is  no  great  concern  over  how  large  the  matrices  may  be  because 
of  the  capability  of  modern  computing  machines.  Matrices  as 
large  as  10  x  10  or  10  x  20  may  present  no  difficulty,  for  example; 
this  point  has  a  bearing  on  whether  average  values  of  w  within  an 
interval  are  used  or  not.  Thus,  if  large  matrices  can  be  handled 
without  difficulty,  it  may  be  as  good  or  better  to  use  a  large 
matrix  and  no  averaging  of  w  as  contrasted' to  the  use  of  small 
matrices  and  average  w  values.  A  point  of  uncertainty  on 
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matrix  size  should  be  mentioned,  however.  Because  experience 
with  the  schemes  has  not  yet  been  obtained  on  a  broad  basis,  we 
don't  Know  whether  difficulties  with  numerics  will  be  encountered. 
We  don't  know,  for  example,  whether  the  large  matrices  may  be 
ill-conditioned,  thereby  making  inversion  difficult.  The  evidence 
provided  by  the  example  cases  of  this  paper  indicates  the 
matrices  are  not  ill-conditioned,  but  whether  this  is  also  so  for 
very  large  matrices  has  to  be  established.  In  spite  of  these 
questions  on  inversion,  It  is  to  be  noted  that  the  developments 
given  in  the  preceding  section  show  £hat  aeroelastlc  problems  may 
be  solved  without  inversion  of  the  K  matrix;  thus,  the  question 
of  matrix  inversion  may,  in  reality,  not  be  of  concern. 

We  conclude  by  identifying  in  summary  the  pertinent  equations 
for  w  which  require  numerical  integration.  These  are  equations 
(38),  giving  w  for  a  subsonic  dipole,  and  equations  (44),  (57), 
and  (62)  which  represent  the  kernel  function.  These  equations, 
expressed  in  a  convenient,  nondimens lonal  form,  should  be  evaluated 
for  various  K  =  <oc/2U  and  M  values.  A  tabulation  of  results 
would  provide  a  universal  set  of  numbers  and  would  obviate  the 
necessity  of  evaluating  the  expressions  in  each  new  problem 
treated.  Application  of  the  techniques  advanced  in  this  paper 
to  solve  various  oscillating  flow  problems  would  be  a  rather  simple 
routine  task. 
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,  as  prescribed 


Figure  1 


Regions  of  Concern  and  Boundary  Conditions 


Concentrated  load 
'  Downwash  point 


Grid  Pattern  Locating  Concentrated  Loads  and 
Downwash  Points 


Slice  a 


"averaged" 


Slice  b, 
"averaged " 


(b)  Averaged  velocities 


Figure  3.  Vertical  Velocity  Field  for  a  Narrow  Horseshoe  Vortex 
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Figure  4.  Basic  Notions  of  Scheme  II 
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Figure  8.  Downwash  Due  to  a  Unit  Load  as  Obtained  from  Wake  Dipole 


Figure  9*  Downwash  Values  for  Solution  3 


Figure  10.  Downwash  Solution  for  Unit  Load  by  Finite 
Difference  Technique  with  z  Finite 


Figure  11. 


Infinite  Array  of  Equally  Spaced  Concentrated  Loads 
Replacing  a  Continuous  Line  Source 


Figure  12.  Example  Airfoil  Cases  Treated 
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(c)  Equal  spanwlse  intervals 


(d)  Spanwlse  Intervals  vary  proportional 
to  chord  intervals 
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